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Top quark pair production close to threshold at a next linear col- 
lider is discussed, with special focus on the influence of rescattering 
corrections. Numerical results are shown for the differential cross- 
Q -1 ' section, i.e. the momentum distribution and the forward backward 

c| , asymmetry, and for moments of the angular distribution of leptons 

arising from the semileptonic top quark decay. 

X 

1 Introduction 

The top quark certainly is one of the most interesting of the presently known 
elementary particles, partly because it is the one we know least of, but also 
because its large mass implies that we can use it for a variety of unprecedented 
QCD studies. A future 500 GeVe + e~ collider will provide a clean laboratory 
where we can investigate the top quark's properties to a high precision by 
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performing a detailed analysis of its production close to threshold. It is thus 
not surprising to find a long series of papers concerned with this subject, i.e. 
leading to the following conclusions: 

1. The excitation curve cr(E), i.e. the total cross-section for it-production, 
allows for a precise measurement of a s and the top quark mass. This 
is in analogy to the charmonium and bottomonium systems with the 
difference that perturbative QCD is far more reliable in the ti system 
because of the considerably higher scale and, as will be explained below, 
the large top quark width. 

2. The momentum distribution da/d\p\ and the forward-backward asym- 
metry will provide a second, independent determination of a s and also 
measure the top quark width T t and consequently give us a handle on 
V t h or physics beyond the standard model. 

3. The size of the top quark polarization, which is predicted to amount 
to ~ —40% even for unpolarized beams, is governed by the top quark's 
couplings to the Z-boson and the photon. Hence its measurement will 
give us information about these presently unknown parameters. 

We will also gain some information about the QCD potential, which deter- 
mines the dynamics of the nonrelativistic "toponium" system. 

When compared with charmonium and bottomonium there seems to be 
a big problem connected to the top quark, caused by its enormous width of 
about 1.5 GeV: the heavy quarks will decay before resonances can build up 
and thus we won't see narrow peaks in ct(E). Furthermore, the width cannot 
be neglected in calculations. 

But as Fadin and Khoze pointed out several years ago [0 , this obstacle in 
fact turns out to be an advantage, because it implies that top decays before 
hadronization effects can spoil the purely perturbative behaviour. The width 
acts as an infrared cutoff and protects us from the poorly known physics at 
large distances. In addition, the study of momentum distributions and top 
polarization is only made possible by the dominance of the single quark 
decay. J/ ip and T for example mainly decay via qq annihilation and thus all 
information about momentum distributions is lost. 

In the following section I will briefly describe the calculation of the cross- 
section for tt production including terms of the order (3. These results are 
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however inconsistent without the inclusion of 0(a s ) corrections, which I will 
present in section 3 for the differential cross-section and in section 4 for the 
top quark polarization. I will not dwell too much on the technical parts of 
the calculation. The reader who is interested in more details is advised to 
read @. 

2 "Born-level" results for tt production 

Close to threshold it is convenient to employ a nonrelativistic approximation 
and use the corresponding variables: the top quark three-momentum in the 
laboratory frame p, with p = |p| and (3 = p/m t , and the energy measured 
from nominal threshold E = ^fs — 2m t = 0(f3 2 ). We will also need the angle 
d between the top and electron beam directions and the three unit vectors 
n., , n ± and n N which are defined in Fig. [I]. 




Figure 1: Orientation of the reference frame. 



In the kinematic regime j3 <C 1 one usually encounters the problem that 
Coulomb gluon exchange is not suppressed because a B / (3 is not a small pa- 
rameter and thus all terms of the form (a s /f3) n have to be summed up. This 
is usually solved by using the Schrodinger equation to calculate the wave- 
functions for the bound states and expressing the observables in terms of the 
latter. In the ti system there is however the additional complication already 
mentioned: the large width which cannot be neglected. As a consequence 
there will be no isolated resonances in the cross-section and thus the usage 
of wave- functions is very inconvenient, because lots of them would have to 
be summed, including interferences. 

The idea of Fadin and Khoze IJ to circumvent this was to use the Green 
function instead, which by construction is a sum over all wave-functions: 



G(r,r^)= £ 



.i*. n E - E n + ie 
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Figure 2: Graphical representation of the Lippmann-Schwinger equations. 

A method to see this function appear in a diagrammatical approach in 
momentum space is as follows]]: one can express the fact that all gluon 
exchanges have to be summed (in ladder approximation) in the form of self- 
consistency equations for the vertex functions Tq, with C = 7^ or -y^ 5 in the 
case of interest for e + e~ annihilation, 

r c = C + J-0^ 4 (4 4 ™s) D^ip-k)rS F (k+^)r c (k,q)S F (k-^, (1) 

which are also shown in Fig. |2|. In the nonrelativistic approximation, ne- 
glecting terms of the order /3 2 , they can essentially be reduced to the integral 
equations^] 

^ , ^ f d 3 k P'k Tr/ , , /Cp(k, E) , , 

/C P (p,E) = 1+ / Tn ~ 2 V{p-k) - P l\ J - . (3) 



(2vr)3 * J E-^- + i T 



mi 

The functions K, in turn are tightly connected to the Green function Q which 
solves the Lippmann-Schwinger equation 

2 - d 3 k 



E-^- + iT t 

m t 



G(p, x, E) = + J j^y 3 V(p - k)G(k, x, E) : (4) 



expanding Q in a Taylor series around x = we find 

£(p, x, E) = G(p, E) + x • ip F(p, E) + ... 



2 This method has hrst been used in 
3 See [pi for details of the derivation 
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with 



G(p, E) = G (p, E)JC s (p, E), F(p, £) = G„(p, E)K P (p, E) (5) 

and 

G (p, E) = 1 . (6) 

m t 

G and F thus are the S- and P-wave Green functions, respectively. They form 
the nontrivial dynamical ingredients for the production cross-section which 
is otherwise completely determined by electroweak coupling constants. For a 
realistic potential they have to be determined numerically, and it turns out 
convenient to perform this calculation directly in momentum space using a 
method introduced in ||. 

We will employ the following conventions for the fermion couplings 

v f = 21} - 4q f sin 2 9 W , a f = 21} (7) 

and the useful abbreviations 

«i = (QeQt + v e v t d z ) 2 + (a e v t d z ) 2 a 2 = 2a e v t d z (q e Qt + v e v t d z ) 

a 3 = 2a e a t d z (q e q t +2v e v t d z ) a 4 = 2a t d z (q e q t v e + (v 2 + a 2 e )v t d z ) (8) 

J _ s 1 

Z {s~m 2 z ) 16 sin^ B w cos 2 W ' 

P± denotes the longitudinal positron/electron polarization and the variable 
X is defined as 

With these parameters the differential cross-section for the production of a 
ti pair with a top quark of three- momentum p and spin projection 1/2 on 
direction n can be written in the form 

d 3 a 1 da / . . > \ 1 + P • n , , 

^ = w*( 1 + 2 ^ (p ' £)co ^)^— (10) 

T = ^{l-P + P-){*i + Xa2)-T t -\p-G(p,E)\ 2 . (11) 

^4 FB denotes the (momentum dependent) forward-backward asymmetry 
and P can be interpreted as top quark polarization vector. These quantities 
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are given by 



A FB ( V ,E) = -CM? R (v,E) (12) 

P,(p,E) = Cj + Cj.^p^costf (13) 

P ± (p,E) = C ± -^ R (p,E)sin^ (14) 

P N (p,£) = C N -^( P ,E)sm^ (15) 



with P ||: _l iN = P • n ||j±)N , 

C»(x) = (16) 

°,W = ' 1 -^ (ai + ^ < 17 > 

Cx(x) = (18) 

c » w = SSTT^" (19) 

and (p R = Keip, tp l = Imip. They have first been derived for polarized beams 
in ||, apart from P N which has been supplemented in |J. The forward- 
backward asymmetry has been calculated first, for unpolarized beams, in 

1- 

The coefficient functions C x are displayed in figure |^. It is not difficult 
to see that for small f3 the top quark spin is mainly aligned with the elec- 
tron beam direction, apart from corrections of order (3 that depend on the 
production angle 

However, the size of these O{0) corrections is not really given by the 
expressions above, although they are correct in a formal sense. Close to 
threshold (3 ~ a s and consequently 0(a s ) terms might be as large, or put 
another way, they should be included in the 0((3) analysis. 

Some of these QCD corrections are not too difficult to implement because 
they are universal. The jit and Zti vertex corrections caused by the exchange 
of transverse gluons for example can be incorporated to the required order 
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Figure 3: The coefficient functions C x for \/s/2 = m t = 180 GeV. 



by simply multiplying with the following factors |T0| : 



G(p,E) -> (i-^GCp,^ (21) 
F(p,£) — (l-^)F(p,F) (22) 

resulting in a corresponding modification in the definition of ip(jp,E). Cor- 
rections to the individual top decays can be accounted for by using the 1-loop 



result for I\ [11 



There is however a type of corrections which is more difficult because it 
connects the different stages of the production and subsequent decay process. 
I will call them "rescattering corrections", the corresponding diagrams are 
depicted in figure |] (in principle there is a third diagram with gluon exchange 
between the two 6-quarks, which can however be neglected to the required 
order 0). 

Fortunately, their contribution to the total cross-section vanishes || PH , 
but as has first been shown in 0], the contribution to the forward-backward 
asymmetry does not (although it was found to be small in a numerical study 
[0]). And there is no a priori argument why they should be irrelevant for the 
polarization vector, especially for those components which start at order (3, 
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a) 



b) 



Figure 4: Rescattering diagrams. 



i.e Pj_ and P N . For this reason I will discuss rescattering corrections in the 
following two sections. 

3 Corrections to the differential cross-section 

By making use of the nonrelativistic approximation, the calculation of the 
rescattering corrections to the differential cross-section can essentially be 
performed along the same line as that of the leading order result. Therefore 
I will skip that part and only discuss the result. The reader who is interested 
in the technical aspects should consult or ||. 
Defining the two convolution functions 




(24) 
(25) 



(23) 



the corrected formula for d 3 a/dp 3 reads 




(26) 
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|pG(p)| a (l+i|>i) |pG(p)| 2 (l+ ¥l ) 




p[GeV] p[GeV] 

a) E=-3GeV b) E=5GeV 



Figure 5: Shape of the momentum distribution of the top quarks for 
E = 3GeV(the position of the would-be IS resonance) and E = 5GeV. Dashed 
line: no rescattering corrections included; solid line: rescattering contribu- 
tion using the full potential included; dotted line: rescattering contribution 
included using a pure Coulomb potential with a s = 0.187. 



with 

42 (P, E) = Ab(p, E) + C° • ^ R (p, E) (27) 

and y = m^/m^. Rescattering thus has two effects: 

1. There is a shift in the momentum distribution, as visible in the factor 
l+ , 0i(p, E), which has a simple physical explanation: if t decays before 
t, its decay product b attracts the top quark and in this way reduces 
the measured momentum. Figure [5| displays the size of this effect. The 
dashed lines show the leading order result, the solid and (hardly visible) 
dotted lines are the corrected curves, with the correction implemented 
in two different ways: the solid line was calculated using the same po- 
tentialQ in the determination of the Green function G and the function 
ipi ; the dotted line was calculated using a pure Coulomb potential in Eq. 
[23| , with the coupling constant fixed to ct s = 0.187 in order to reproduce 
the solid line. This was done to compare with || where a pure Coulomb 
potential was used as well, but with a s = a s (a s (M z )m t ) ~ 0.15. Of 

4 Its actual expression can be found in || 
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10 20 30 40 50 60 70 80 10 20 30 40 50 60 70 80 

p[GeV] p[GeV] 

a) E=-3GeV b) E=5GeV 



Figure 6: Rescattering corrections to the forward-backward asymmetry. 
Dashed line: S-P-wave interference contribution; solid line: full result. The 
little arrows indicate the position of the peak in the momentum distribution. 



course the scale governing rescattering corrections is of the order of 
the top momentum p ~ a s (Mz)m t , but when using a fixed coupling 
constant, its actual value is nevertheless unknown a priori. 

2. There is a change in the forward-backward asymmetry which depends 
on the beam polarization, as can be seen from Fig. [f| The dashed 
line again shows the uncorrected, the solid line the corrected result. 
The small arrows indicate the position of the peak in the momentum 
distribution, i.e. the region which is most relevant. The reason why 
the corrections were found to be small in || is that these authors only 
considered unpolarized beams. 

There are also a few things worth noticing: 
The role of T t as an infrared cut-off is obvious in Eq. |24|. For vanishing width 
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the integrand would behave as 



nik- P i) i 

|k — p| |k — p| 3 

for k ~ p, but the factor in brackets weakens this behaviour considerably: 

V(\k-p\) f Ft |k-pk 1 1 

1 ~ tt: — — arctan — ^ — ) ~ r> |k-p|<r t . 



|k- P | v | k _ p | r t > 3r t 2 k 

The corrections to the total cross-section indeed vanish because they are 
proportional to the integral 



/^|G(p,£)| 2 ^(p,£) = 0. 



This can be seen by inserting Eq. |23] and interchanging the integration vari- 
ables p and k. 

It is advantageous to solve the Lippmann-Schwinger equations in momen- 
tum space because this allows for a relatively easy numerical determination 
of i/ji and ip2- I n the coordinate space approach additional Fourier transfor- 
mations would be necessary leading to a loss in accuracy. 



4 Top quark polarization 

The calculation of the rescattering corrections to top quark polarization is 
more involved than that of the spin-averaged differential cross-section. In the 
Born-level calulation the quark is essentially treated as a stable particle, the 
free Green function Go being the only place where its decay enters explicitly 
through the usage of the complex energy E + iT t . The top quark spin is con- 
sidered a well defined quantity and thus the appropriate projection operators 



are used in the derivation of Eq. 10 



This is however no longer possible in the evaluation of the diagram shown 
in Fig. f|(b). There the top quark appears as a virtual particle only and 
thus the notion of its polarization is highly questionable. To find a way out 
one may ask: How will the top quark polarization be measured in a real 
experiment ? 
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The answer to this question is known: the experimentalists will use the 
angular distribution of the leptons from the decay t — > blu, because for a free 
top quark in its rest frame, this is given by 



dT t _ blu {P) = 1 
d cos 9 2 



r t _> M „(l + |P|cos0 



where P denotes the top quark polarization vector and 9 the angle between 
P and the lepton momentum. 

The generalization of this formula to tt production in e + e~ collisions 
would be 

/F>rr A'irr 1 /TV. ,, (T>\ 

(28) 



d 6 a d 3 a 1 dIW,(P) 



dp 3 dEi dQi dp 3 T t dE l dtt t 

which expresses the fact that first a top quark with momentum p and polar- 
ization P has to be produced, which may then decay semileptonically. Di- 
agram ||(b) however destroys this factorization because the production and 
decay stages are no longer separated. 

Of course the left hand side of Eq. ^ could be calculated nevertheless, 
at least in principle, however it would lead to a rather complicated result. 
A method to obtain simpler expressions involving an easier calculation is to 
consider moments of the lepton spectrum as defined by 

. /d 3 a\-i f ,^ da(e + e~ — > blvbW") . 

< n ^M J dE > da ' wd^dE, <*> < 29) 

or 

= f***kM, dpdnpdn , dEl - <"'), (30) 

where / is the lepton four-momentum and n may be any convenient four- 
vector, for example (in the ti rest frame) 

9T, H = (0, n N ) n ± = (0,n ± ) (31) 
n N = (0,n N n (0) = (1,0,0,0). (32) 

This method has several advantages [[|: 

• it always leads to well-defined, physical, measurable quantities, 
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it is far easier than calculating the full multi-differential cross-section 
because the leptonic phase space can be integrated out completely, even 
in one of the first steps, 

it still gives us the same information on P, or more exactly, it can be 
used to define a quantity which is equivalent to P in the absence of 
rescattering. Neglecting QCD corrections one easily derives from (f28|) 

1 + 2y + 3y 2 



(nl) BoIn = BR(t - Uv) | ( | ^ 

where t and P are four-vectors given by 
t = (m t ,p) . P = ( 



[n-t) + ^(n.P) 



— C°cos^,P) 



(33) 



(34) 



The appearance of the 0-component of P is welcome, because if P is to 
be interpreted as a polarization vector, it has to fulfill the transversality 
condition t ■ P = 0, which indeed it does. Remember that we are work- 
ing in the ti rest frame instead of the t rest frame (the corresponding 
change in P is of the order {3 2 and thus does not appear here). 

Rescattering modifies equation fl33|), the corrected expression reads 

1 + 2y + 3y 2 



(nl) 



BR(t Uv) 



4(1 + 2y) 



mt, 



(n ■ t)(l + e) + -^(n ■ P + n ■ 5P) 

o 



with 



v0--y) 



l + 2y)(l + 2y 
2 + 3y-5y 2 - 



-3y 2 ) 
12y 3 



C> R (p, E) costf 



(35) 



(36) 



l-2y 



0\2 



'l + 2y)(l + 2y + 3y 2 
1 - Ay + 3y 2 



5P 



5P, 



4(1 + 2y + 3y 2 
3(1-3^ 



1 



l + 2y 

3 cos 2 tf)C°V> 3 (p,£) 



^ R (p, -E") cos-# 



(37) 



2(1 + 2y + 3?/ 2 , 
3(l-4y + 3y 2 

8(l + 2y + 3y 
0, <5P = 0. 



^ R (p, E) sin^ 
C> 3 (p,^) sin 2d 



(38) 
(39) 
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Figure 7: Subleading part of the longitudinal component of the polarization 
vector P for i? = 0. The solid line shows the complete result for unpolarized 
beams, the dotted line the S'-P-interference contribution. The dash-dotted 
line shows the complete result for fully polarized beams, in which case the 
S'-P-wave contribution vanishes. The arrows indicate the position of the 
peak in the momentum distribution. 



The new function ip 3 is defined as 

d 3 k 



G{k, E) arctan ^ 



x 



(2tt) 3 
/ p - (k-p) ^2 
^ p|k-p| ) 



1 



|k- P | 



(40) 



and its presence destroys the otherwise nice similarity between P and 5P. 
Numerically its contribution is however negligible for realistic top quark 
masses. The same applies to e so that at least in practice the form of the 
Born-level result is recovered. One should note that there was no a priori 
reason to expect this form to remain the same. 

In addition I should mention that ipi, ip 2 and ip 3 could be interpreted as 
additional S-, P- and P-wave contributions, a viewpoint motiviated by the 
way they are defined and the way they enter the expressions. 
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p[GeV] p[GeV] 

a) E=-3GeV b) E=5GeV 

Figure 8: Transverse component of the polarization vector for d = tt/2. 
The solid and dotted lines correspond to unpolarized beams, the dash-dotted 
and dotted lines to fully polarized beams. The solid and dash-dotted curves 
show the full result, the dashed and dotted the S-P-interference contribution. 
The arrows indicate the position of the peak in the momentum distribution. 

Figures [7| and || show the predictions for the subleading term of the longi- 
tudinal polarization and for the transverse part, respectively. The solid and 
dash-dotted lines correspond to the full result for unpolarized/fully polarized 
beams, the dashed and dotted lines to the 0{(3) term only. It is quite evident 
that rescattering has a strong influence on the result, which is due to the fact 
that the 0(/3)-terms are suppressed by ratios of electroweak couplings ap- 
pearing in the coefficients C* and C±. <£> R and ip n are of comparable size as 
they should be because (3 ~ a s . 



5 Summary and conclusions 

In this paper I have discussed recent results on observables in top quark pair 
production near threshold at e + e~ colliders. I have sketched the calculational 
procedure used to derive the predictions for the differential cross-section and 
top quark polarization for arbitrary beam polarization which are correct up 
to terms of the order p 2 / m t because they also include the rescattering cor- 



15 



rections. The method of moments proved to be a very useful tool, which 
could of course also be used for other production mechanisms or far above 
threshold. 

It turned out that rescattering corrections are not generically small and 
are indeed required for a consistent prediction of the components of the top 
quark polarization vector lying in the production plane. They are of minor 
importance for the forward-backward asymmetry but might nevertheless be 
visible in precision studies, especially if polarized beams are used. An impor- 
tant result was that there are no corrections to the normal polarization. This 
component is of special interest because it is the place where non-standard 
CP violating interactions would manifest themselves. In connection with the 
leading term of P.. , which can be measured by averaging over the top direc- 
tion, i.e. via ({n.. -P)), it provides an opportunity to determine the top quark 
couplings ai . . . a 4 . 

I would like to express my gratitude to the organizers of the XXXVI Cra- 
cow School of Theoretical Physics, especially M. Jezabek, for their invitation 
and for providing such a nice atmosphere. I would also like to thank J.H. 
Kiihn for carefully reading the manuscript and giving useful advice. 
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